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In two intersecting many-particle streams, one can often find the emergence of oscillatory pat- 
terns. Here, we investigate the interaction of pedestrians with vehicles, when they try to cross a 
road. A numerical study of this coupled pedestrian-vehicle delay problem has been presented in a 
previous paper. Here, we focus on the analytical treatment of the problem, which requires to use a 
simplified car-following model. Our analytical results for the phase transition to oscillatory pedes- 
trian and traffic flows and the average waiting times are well supported by numerical evaluations 
and give a detailed picture of the collective dynamics emerging when pedestrians try to cross a road. 
The mathematical expressions allow one to identify the dependence on model parameters such as 
the vehicle or pedestrian arrival rate, and the safety factor of pedestrian gap acceptance. We also 
calculate a formula for the vehicle time gap distribution, which corresponds to the departure time 
distribution of a M/D/l queue. 

PACS numbers: 89.40.+k,47.54.+r,47.55.-t, 02.50.-r 
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I. INTRODUCTION 

Pattern formation is a wide-spread feature of driven 
many-particle systems. In particular, oscillatory patterns 
are found in fluids, granular materials, colloidal systems, 
and traffic flows. One typical example are stop-and-go 
waves in traffic flows on freeways caused by a delayed 
adaptation to changing traffic conditions 0, 0, 0, 0> HJ . 
Emergent oscillations have been discovered in so different 
systems as the density oscillator ticking hour glass Q , 
RNA Polymerase traffic on DNA §| , pedestrians passing 
a bottleneck 0,0] , or ants 0] . Oscillatory patterns have 
also been found in two intersecting pedestrians streams 
[l2| or simulations of colloidal systems [l3|. 

Although the subject is rather old 0, 0, 0, the 
crossing of vehicle streams by pedestrians has recently 
attracted an increasing interest, also among physicists 
[II EH El However, the problem of interactions 
between vehicles and pedestrians, when pedestrians are 
trying to cross a road, has not yet been sufficiently un- 
derstood. The mathematical investigation of this problem 
will be the subject of this paper. Numerical studies have 
shown a transition from crossing the road one by one 
or in small groups to coupled oscillations of pedestrian 
and vehicle flows, if pedestrians use small gaps to cross 
the road [2lJ . In the following, the dynamics of this phe- 
nomenon and the parameter-dependence of the transition 
point will be investigated analytically. 

Our paper is organized as follows: Section ^ formu- 
lates the model for the pedestrian and vehicle behavior 
and their interactions. Moreover, we calculate a formula 
for an idealized vehicle time gap distribution. In Sec. IIIII 
we will derive analytical results on the dynamic behavior 
of interacting pedestrian and vehicle flows. Moreover, we 
will compare these results with numerical evaluations of 
computer simulations of the underlying model. Our ana- 
lytical formulas for the transition point and the waiting 
times of pedestrians and cars are well compatible with 



numerically determined data. Finally, we will summarize 
and discuss our results in Sec. II VI which are relevant for 
many systems with intersecting flows or competing pro- 
cesses. 



II. FORMULATION OF THE MODEL 
A. Vehicle Behavior 

In our simplified model of vehicle dynamics, cars are 
treated as moving objects of length / - We assume a con- 
stant arrival flow Qarr of vehicles and that new cars try 
to enter the investigated road section with a probability 
q = Qarr dt per time step dt. This implies an exponen- 
tial time gap distribution, which is modified by vehicle- 
vehicle interactions (see Sec. Ill B|) . In fact, a vehicle with 
the speed v following a leading vehicle with speed v* is 
assumed to decelerate with dv/dt = —a, if v > and 
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where Ax denotes the distance between the two vehi- 
cles, Iq the vehicle length and do the preferred minimum 
bumper-to-bumper distance among cars. This condition 
guarantees accident-free driving [22. For a ">"-sign in 
Eq. 0J, the vehicle accelerates with dv/dt — a, delayed 
by the reaction time T, until the maximum (free) speed 
vo is reached. For an "="-sign in Eq. (JTJ, the veloc- 
ity is not changed, i.e. dv/dt = 0. The above contin- 
uous car-following model may be called the constant- 
deceleration-delayed-acceleration model (CDDA model) 
and has some similarities with the slow-to-start cellular 
automaton model 23] . 

We assume that pedestrians enter the street at the 
crossing point O, when they consider it safe (see 
Sec. Ill C|) . Moreover, crossing the road takes a time pe- 
riod t. In order to avoid accidents with pedestrians, vehi- 
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cles decelerate with dv/dt = —a if necessary. We consider 
two different deceleration rules: 

a) Careful drivers: The closest car to a pedestrian 
on the street decelerates, if the distance d(t) — 
—x(t) to the crossing point O is within the range 



< d(t) <d + 



2n 



(2) 



where do is the safety distance that a car should 
keep from a crossing pedestrian. We assumed this 
safety distance to be identical to the minimum 
bumper-to-bumper distance among vehicles ap- 
pearing in ifljl. 

b) Aggressive drivers: The closest car starts to de- 
celerate at the time to determined so that the dis- 
tance to the pedestrian corresponds to the safety 
distance d(t n + r) — do at the time t n + r when the 
last (the n-th) pedestrian on the street (entering at 
time t n ) leaves the road after the crossing time r. 

After the last pedestrian has left the street, i.e. at time 
t n + r, the car accelerates with dv/dt = a, until it has 
reached its desired velocity v again. The characteris- 
tic distance between stopped vehicles in a queue is the 
vehicle length lo plus the minimum bumper-to-bumper 
distance do, which defines the jam density 
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In the following, we will assume that a car starts to ac- 
celerate after its leader delayed by the reaction time T. 
This implies that the following car has reached the posi- 
tion of the leading car in the queue after a time period 
T + \/2(d + lo)/a and that the distance to the leading 
car is lo + do + Tvo, when the following car has reached its 
maximum velocity Vq. Therefore, the outflow from a traf- 
fic jam starts with a value of (T + ^2(d + ^o)/o) _1 an d 
eventually reaches the characteristic (maximum) value 
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while the traffic jam (queue) resolves upstream with the 
characteristic speed 
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due to the distance lo + do between queued cars and 
the delay T in acceleration. Moreover, when a vehicle 
is stopped at point x(t) = —do, the forming traffic jam 
behind it propagates upstream with the velocity |24j 
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which depends on the vehicle arrival rate Q 3 



The proposed simple car-following model essentially 
reflects the features of the section-based, fluid-dynamic 
traffic flow model proposed in Ref. |24|], with the only 
difference that the acceleration and braking processes re- 
quire time periods of T + Vo/a and T/vo, respectively. 
Apart from scaling time and space variables in order to 
get rid of two more model parameters, it is hard to think 
of any further simplification of the above vehicle model 
without sacrificing fundamental properties of traffic flows 
such as the constant outflow from traffic jams and the 
characteristic jam resolution speed |2j|. Nevertheless, it 
may be interesting to study the limit a — » oo of unlim- 
ited acceleration possibilities, which eliminates accelera- 
tion and deceleration times. More realistic variants of the 
above car-following model, however, should distinguish 
different acceleration and deceleration strengths a and b, 
which have been set equal here for the sake of simplicity. 
A stochastic variant of this model describing a fluctu- 
ating acceleration behavior would be also interesting to 
study. 



B. Idealized Vehicle Distance Distribution 

In our vehicle simulations, we have generated vehicles 
with initial velocity v — at the upstream boundary of 
the simulation stretch according to the exponential time 
gap distribution Qarre - ^"" T , where T" denotes the ac- 
tual time gap. However, according to our car-following 
model, vehicles had gained at least their preferred dis- 
tance D = lo + do + VqT, when they reached the maxi- 
mum speed vq- According to theoretical considerations, 
this changed the effective time-gap distribution at the 
crossing point to 

P(T') = Q arr T S(T' - To) 

+ (l-Q a rrT )Q a „e- Q ^ T '- T ^e(T' -T ) (7) 

with To = D/vo (see Appendix), when no vehicles at the 
entry point were dropped. That is, a fraction Q air To of 
vehicles will follow with the desired time gap To, while 
the rest has an exponentially distributed, larger time gap 
T' > Tq. S(y) denotes Dirac's delta function, while the 
Heaviside function Q(y) is 1 for y > and otherwise. 

Our exponentially distributed vehicle generation mech- 
anism sometimes causes a virtual queue of vehicles at the 
upstream boundary, which can be avoided by generating 
vehicles according to the resulting time gap distribution 
(|7|). In fact, our implementation of the boundar y co ndi- 
tions corresponds to a M/D/l queuing system |26l l27]|. 
i.e. to a queue with Poissonian distributed Markovian 
arrivals (where the time gaps between successive arrivals 
are exponentially distributed), while the service rate 1/To 
is assumed to be deterministic. (The "1" stands for one 
"channel", i.e. no parallel service.) 

Now, let Po be the probability that no vehicle is wait- 
ing in the queue to be served, i.e. to enter the road. The 
probability of releasing the next vehicle with a time gap 
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T' = T is then given by the probability (1 — P ) of 
having queued vehicles waiting to enter, plus the prob- 
ability Po(l — e _( 5 a,rTo ) that we have the no-queue case 
and a vehicle arrives during the service time To. In cases 
with no queue where the time gap T' of the next arriving 
vehicle is greater than To, we have an exponential time 
gap distribution Q arr e~® a " T /e~^ arrTo , where e -< 2 arrT ° is 
the normalization factor of the conditional probability of 
finding time gaps larger than To. Altogether, we obtain 
the time gap distribution 



P(T') = [(l-T )+To(l-e- Qa " To )]<5(T'-T ) 
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i.e. that the vehicle flow Q arr and, therefore, the average 
time gap remains unchanged, we find 



TO = (1 - QarrTo) e Q *" T ° . 



(10) 



This implies the idealized vehicle time gap distribution 
0, which will be necessary to evaluate the expected 
waiting time of pedestrians for a suitable time gap to 
cross the road (see Appendix). 



2) Daring pedestrians enter the road, if a car with 
velocity vp would not come closer than the preferred 
safety distance dp, if it decelerated with dv/dt = 
—a in order to avoid an accident. This implies the 
reduced safety factor 
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In this case, a single pedestrian can force a car to 
stop, namely when entering at a vehicle distance 
d(t) = d +v 2 /(2a). 

Realistic values of the safety factor a are expected to be 
above <J2- 

For the following analysis, we will identify the time 
point t = with the time when the first pedestrian(s), 
who cause(s) a vehicle to decelerate, enter(s) the road. 
The entering time of the next entering pedestrian is de- 
noted by t\ , the entering time of the /cth following pedes- 
trian by ifc, and the entering time of the last (nth) fol- 
lowing pedestrian before the car passes point O by t n . 



D. Simulated Dynamic Behavior of Interacting 
Vehicle and Pedestrian Flows 



C. Pedestrian Behavior 

We will assume that pedestrians enter the sidewalk of 
the street at the crossing point O with probability p = 
A dt per time step dt, i.e. A denotes the arrival rate of 
pedestrians. If there is no sufficient gap in the vehicle 
stream to cross, they accumulate around point O, but 
they start immediately to enter the road at time t, if 
v(t) = (i.e. if the vehicle velocity is zero) or if 



d(t) > d and At(t) := 



d(t) 

v(t) ~ 



(11) 



(i.e. if the distance d(t) is larger than the preferred safety 
distance do and the time gap At is large enough to cross 
the road). Here, At is the time to collision of the nearest 
approaching vehicle and a a safety factor of pedestrians, 
r is the time period required for a pedestrian to cross 
(one lane of) the road. We may distinguish two limiting 
cases of gap selection, i.e. interactions with approaching 
vehicles: 

1) Careful pedestrians assume that cars may not 
decelerate and approach with their desired velocity 
vp. They cross the road only, if the car at no time 
comes closer than the preferred safety distance do, 
which implies the following choice of the safety fac- 
tor: 

dp 

v t ' 
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Simulations of vehicles interacting with pedestrians 
crossing a street have recently shown an interesting phe- 
nomenon: While for large enough values of the safety 
factor a, pedestrians cross the road one by one or in 
small groups, one finds alternating pedestrian and ve- 
hicle streams if the safety factor is smaller than some 
critical value ap- This value can be exactly calculated for 
the above model (see Eq. (|33fl ). which shows qualitatively 
the same dynamic behavior like the variant of the IDM 
model studied in a previous publication |2lT | . Representa- 
tive simulation results for the above proposed pedestrian 
and vehicle model are displayed in Fig. ^ The parame- 
ter values used in this paper are a — 1 m/s 2 , r = 2 s, 
T = 0.9 s, /o — 4 m, do — 2 m, and vp = 15 m/s, and our 
numerical investigation focusses on careful drivers. 

The reason for the observed oscillations is that pedes- 
trians can force vehicles to stop, if they choose small time 
gaps At. However, if vehicles are stopped, they have to 
wait until there is a gap of period r or larger in the 
pedestrian stream, before they can accelerate again. Dur- 
ing this waiting time, a vehicle queue is formed, which 
can become very long, dependent on the vehicle arrival 
rate. Pedestrians cannot cross the road again, before this 
queue is completely dissolved, at least if 
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i.e. if the time gap between successive vehicles having left 
the queue is too short for pedestrians to enter the street, 
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and if 
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FIG. 1: (a) Representative space-over-time plot of vehicle tra- 
jectories for careful drivers and the pedestrian safety factor 
a — 1.05. Pedestrians may stop cars, which causes vehicle 
queues. These suppress the crossing of newly arriving pedes- 
trians until the vehicle queue has completely dissolved, (b) 
Representative space-over-time plot of vehicle trajectories for 
the larger safety factor a = 1.25, for which pedestrians use 
large gaps only. As a consequence, pedestrians do not stop 
cars completely when they cross the street, and no vehicle 
queues are formed. 
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i.e. if the time gap with respect to the second car in the 
queue at the time \/2(l + dp)/a > T (when the back of 
the first vehicle has passed the crossing point O) is not 
large enough for pedestrians to enter the street. 

In summary, we may have alternating time periods in 
which pedestrians can cross the road and time periods in 
which cars can pass point O. In the following sections, 
based on statistical approaches, we will try to estimate 
the time period until a sufficiently large gap in the ve- 
hicle flow occurs to allow pedestrians a crossing of the 
road. Likewise, we will calculate the time period until 
queued vehicles find a large enough gap between cross- 
ing pedestrians, allowing them to accelerate again. An- 
alytical results can be only gained for simple models as 
the one proposed above. Nevertheless, we expect quali- 
tatively similar relationships for a broad class of other 
traffic models. 



III. ANALYTICAL RESULTS AND 
COMPARISON WITH COMPUTER 
SIMULATIONS 

A. Dynamics of Vehicles Reacting to Pedestrians 

Let to be the time point when the car starts to de- 
celerate as response to a crossing pedestrian. According 
to Sees. Ill A1 and ITTHl we find that the time to collision 
evolves in time according to 

Am = m = d{o)-vot = m_ t if f<to . (16) 

v(t) v v 

For careful drivers, i.e. case a), the start time of deceler- 
ation can be determined as 



tn = 



d(0) - d v 



vo 2a 
This yields the time to collision 

_ d + vo 2 /(2a) - v (t - t ) + a(t - t ) 2 /2 
v - a(t - t ) 

vo_ _ t-t d 

2a 2 vp ~ a(t — to) 

(see Fig. EJ and the vehicle velocity 
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if t > t (18) 



v(t + r) =v -a(T-tp) = ^--ar + a ^ d ° (19) 

2 vp 

after the pedestrian has crossed the road. If the vehi- 
cle velocity at the beginning of the braking maneuver 
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FIG. 2: Time- dependent time to collision At(t) = d(t)/v(t) 
for careful drivers [see Eq. 1181 1. when pedestrians would en- 
ter the road with probability p = 1 and a = 1.05 (symbols 
= numerically determined values, solid line = analytical for- 
mulas). Due to the braking maneuver, the time to collision 
goes down in the beginning, but it grows again later on, as 
the vehicle comes to rest at the finite distance d(t) — do to 
the pedestrian. 



is v(to) < vq, one just has to replace vq by v(to). For 
aggressive drivers, i.e. case b), we find 



t(] = T - 



2v t d(0) - d 



(20) 



At(t) = 



and 
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if t > to, (21) 



v(*o +t) = v - a(r - t ) 
= v 



y / 2avoT - 2a[d(0) - d ] ■ (22) 



That is, the greater the initial distance, the later will the 
vehicle start to decelerate and the larger will the resulting 
velocity be. Note that, according to the gap acceptance 
rules of pedestrians outlined in Sec. Ill CI the shortest 
distance to a moving vehicle at which pedestrians enter 
the road, is given by arv. 



free velocity w can be calculated as the distance 2(vq — 
Wmin) 2 / (2a) travelled less, divided by the desired velocity 
vq, which results in 



A<br = 
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If the vehicle is stopped, the time lost by the accelera- 
tion and deceleration process amounts to vo/a. On top 
of this, we have to add the average waiting time t w . This 
can be obtained as follows: If Ati denotes the waiting 
time of the first stopped vehicle, the number of vehicles 
queuing up behind it until the first car in the queue starts 
to accelerate is given by Pj am C Ati. The delay of the last 
vehicle in the queue is the queue length I = C Ati, di- 
vided by the queue resolution speed c. As the waiting 
time between the first and the last vehicle in the queue 
progresses approximately linearly, their cumulative wait- 
ing time is given by 
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(24) 

Moreover, upto the time point when the queue formed 
within the stopping time Ati has resolved, another 
PjamlC/ (c — C) vehicles have joined the queue (cf. For- 
mula (1.48) in Ref. 0). While the waiting time of 
the first of these additional vehicles is approximately 
l/c — C Ati/c (as the one of the last vehicle in the first 
part of the queue), the waiting time of the last vehicle is 
basically zero, which implies a cumulative waiting time 
of 



p jam CAti f C CMi 
:-C c 
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p iam C(Ati) 



c 2 



2 c? -cC 
(25) 

Adding this to (|24|l gives the cumulative waiting time 



«am (Ati) 2 CC 
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(26) 



which grows quadratically in Ati (see Fig. [31). 

Finally, dividing this result by the total number 
CAti[l + C/(c — C)} of vehicles yields a very simple 
relationship for the average waiting time, which is just 
given as the average waiting time of the first and the last 
queued vehicle: 



tw — 



Ati 



(27) 



However, the estimation of the waiting time Ati of the 
first stopped vehicle is rather difficult (see Sec. IIIIE)l . 



B. Average Delay to Vehicles 

Let us denote by i> m in the minimum velocity before the 
car accelerates again. If only one pedestrian obstructs 
the car, we have u m j n = u(r), as calculated above. The 
time delay to the car compared to a movement with the 



C. Determination of the Transition Point to 
Alternating Flows 

The long vehicle and pedestrian queues required for 
pronounced oscillations in the pedestrian and vehicle 
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FIG. 3: Average of the cumulative waiting times t c of vehicles 
as a function of the time period Ati the first vehicle in the 
queue has to wait, for different values of the vehicle arrival rate 
Q arr , see Eq. © (symbols = numerically determined values, 
parabolic curves = analytical formula). 



flows can only occur, if vehicles can be completely 
stopped by pedestrians. This cannot happen, if the safety 
factor a of pedestrians is large enough. For small values of 
a, however, there exists a time point after which the 
safety criterion prohibits a further entering of pedes- 
trians into the road. This time point is given by the earlier 
time fulfilling the critical safety condition Ai(fep) = err. 
Together with the expressions for the times to collision 
in Sec. Illl Al this eventually implies 
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for careful drivers. t+ is the first time point at which 
pedestrians may re-enter the road again, as the time to 
collision At(t) increases close to the crossing point [see 
formula (|18() ]. The car reaches its minimum possible ve- 
locity a time period r after t_, i.e. after the latest enter- 
ing pedestrian has left the road at time t_ + r. With 128|) 
this implies 



v(t— + r) = ar(cr — 1) + y/ (aar) 2 — 2ado 



(29) 



for careful drivers. For aggressive drivers, we have to re- 
place do by d(to) — vq 2 /{2a). To exclude stopped vehicles, 
on the one hand, this minimum velocity should be posi- 
tive, i.e. 
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ar 2 > do . 
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On the other hand, vehicles could also be stopped by new 
pedestrians entering the road at a time t > t+ that lies 
before the time t_ + r at which the last pedestrian has 
left the road. Therefore, in order to avoid the stopping 



of vehicles by multiple crossing pedestrians, we have to 
demand 



t+-*- = 2^(or)* 
which results in 
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Together with condition i|30[l we find that a careful driver 
cannot be stopped completely under the condition 



a > do = max 
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(33) 



At the value a — <7n, we expect a transition from con- 
tinuous pedestrian and vehicle flows to alternating flows 
(see Fig. 0} . 



D. Calculation of Earlier Acceleration 

Due to the statistical arrival of pedestrians with a rate 
A = p/dt, it is likely that the time point t n < t_ of the 
last (nth) pedestrian entering the road is smaller than 
the latest possible entering time i_. We are, therefore, 
interested in calculating the mean value (t_ — t n ) = £_ — 
(t n ) of the time gap t_ — t n , where n is an arbitrary 
integer number. For this, let K = t—jdt be the number 
of time steps between the first entering pedestrian and 
t_ . As the probability that no pedestrian enters in a time 
step is given by r = (1— p), (l—p) K is the probability that 
nobody enters between t = and t = t_, andp(l— p) K ~ k 
the probability that the last pedestrian enters at time 
t_ — (K — k)dt = k dt. The expected value of i_ — t n is 
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(see Fig.|SJ. Therefore, if a vehicle is not stopped, instead 
of at time t_ +r, on average it already starts to accelerate 
at the earlier time 



(<„} +T = £_ + T- (1 - p)dt 
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where the last step of this calculation is based on Eq. lj3*T| 
below. With this result, we can now estimate the ex- 
pected value (f m in) of the minimum vehicle velocity t> m i n 
entering Eq. (J2HJ): 

(v mi n)-v(t-+T) =a(t--(tn)) = a(l-e- xt -)/X. (36) 
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FIG. 5: Average time span t_ — (t n ) between the latest possible 
entering of the street by a pedestrian and the time point when 
the last pedestrian actually enters the street as a function of 
the scaled pedestrian arrival rate At, see Eq. 1361 (symbols 
= numerically determined values, solid line = analytical for- 
mula) . 



The higher velocity compared to + r) given by 

Eq. I|29[) originates from the earlier car acceleration, i.e. 
the shorter deceleration time. 
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FIG. 4: Transition point ao to alternating vehicle and pedes- 
trian flows as a function of the dimensionless parameters 
(a) xi = do/(var) obtained for do £ [1.6 m, 2.5 m], (b) 
xi = Vo/(ar) obtained for a 6 [0.5 m/s 2 , 1.2 m/s 2 ], and (c) 
x-a — do/(ar 2 ) obtained for a £ [0.5 m/s , 1.2 m/s ] in com- 
parison with the lower limit 02 of reasonable safety factors 
[see Eq. I13H 1 (symbols = numerically determined values, solid 
lines = analytical formula). Note that the value of ao is con- 
stant for X3 = X1X2 = do/(ar' 2 ) — const. 



E. Estimation of the Waiting Time of the First 
Vehicle 

If a vehicle is stopped by crossing pedestrians after a 
deceleration time i>o/a, it will have to wait until a time 
gap of duration r in the pedestrian flow occurs. A gap of 
length t = N dt or greater occurs with probability 

(1 - p) N = (1 - P y' dt = [(1 - pY'^f = e- A ^ , (37) 

>> .. ' 



i.e. gap sizes are exponentially distributed, as expected. 
Here, we have assumed ln(l — p) w — p, but the required 
small values of p = Xdt can be reached by sufficiently 
small choice of the time steps dt. In fact, in the following 
considerations, we will study the limit dt — > 0. Therefore, 
we have used the value dt = 0.001 s in our computer 
simulations. 

Now, let ki denote the size of the ith gap Tj = U — 
tj_i (i.e. the number of time steps dt with no pedestrian 
arrival). Then, the expected value for the time period 
until a time gap of length r = N dt or greater starts is 
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given by 

oo N 



N 



E E • • • E + 1 + • • • + fc « + !) 

n=0fe 1= fc„=0 

*(!-!>)''!> (1 -)>)'"!>• (1-p)" 



£*•/*!.../ 



dT„ 



x{Ti 
„-A 



Tn)e _ A(Tl +... + T„) e -Ar 



ri=0 i=l 
oo , 

oo 

T E«(!- 



-ATi 



(1 



-Xt\ 



AtN 



3 -At 



1 - 
A 



,-At 



3 -At 



T s— Vs" with s = 1 -e" Ar 



,1 



,A 



r -(l + Ar)] 



___ 

2 



(38) 



That means, the average waiting time for a gap of size 
t or larger starts to grow linearly with the pedestrian 
arrival rate A = p/dt and quadratically with r as long as 
these values are small, but it grows exponentially with 
At, when this value is large. 

Note, however, that the waiting time is reduced by the 
gap between the time M dt := vo/a when the vehicle is 
stopped and the time t n < t_ at which the last pedestrian 
has entered the street before. Analogously to Sec. IIII Dl 
we can calculate the expected value of this time gap as 



(v /a-t n ) = (1 -p) 



I-(I-P) 
p/dt 



M 



I _ g— Auo/a 

A 



(39) 



since we have (1 — p) — > 1 in the limit dt — > 0. As a 
consequence, the expected value (Af i) of the time period 
Aii the first vehicle in the queue has to wait can be 
estimated as 



1 1 _ p-Auo/a 

(AtO = -[e^-(l + Ar)] 



= i(( N "" " - 2- Ar 



(40) 



(see Fig. |SJ). 



F. Average Delay to Pedestrians 

After a time interval Ati, i.e. a time period r after the 
last pedestrian has entered the road, the first vehicle in 
the queue can accelerate again. The time period available 
to pedestrians for crossing the road is Ati + Vo/a, as the 



1000 




FIG. 6: Average waiting time (Ati) of the first vehicle in 
the queue as a function of the scaled pedestrian arrival rate 
Ar (symbols = numerically determined values, solid line = 
analytical formula). 



time period vo/a required to stop the vehicle is usable 
as well. When the vehicle has started to move again, no 
pedestrian will be able to cross the road until the last 
vehicle of the queue has passed point O (at least if err > 
T). This time period can be calculated as p4| 



Ah 



CAtj 



1 + c/v 

c-C 



The expected value of At 2 is 



(At 2 ) = C(At 1 ) 



(41) 



'-^ (42) 



(see Fig. 0|, where we have also taken into account the 
additional amount yj 2do / a required by a vehicle to get 
from x — — d to point O. 

After the last vehicle in the queue has passed, pedes- 
trians have a chance to find a suitable gap of size err or 
larger. A lower bound of the expected waiting time (At 3 ) 
for the occurence of such a gap is calculated in the Ap- 
pendix. In Fig. [SI we compare the resulting expression 



<n> = 



i 



Qa 



r(T,-T ) 



1 — Qarr^O 



- (1 + QarrT*) 



(43) 



corresponding to Eq. (|A7|) with numerical results, where 
T = l/Qout = T+(l +d )/v Q andT» = aT+(l Q +d Q )/v Q . 
This formula gives the expected waiting time (At 3) pro- 
vided that the pedestrian arrives exactly at the time, 
when the last vehicle in the queue passes the crossing 
point. Otherwise, it is an approximation, which neglects 

1. the effect that pedestrians tend to arrive at the 
sidewalk between two vehicles (so that there is an 
incomplete intervehicle time gap, which must be 
added) , 
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400 




<At,> (S) 



FIG. 7: Average time (Atz) needed to dissolve a vehicle queue 
as a function of the average time (Ati) for which the first ve- 
hicle has been waiting, for various values of the vehicle arrival 
rate Q a rr, see formula JSJ (symbols = numerically determined 
values, straight lines = analytical results). 
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FIG. 8: Average waiting time {Ats} until a pedestrian enters 
the road after a vehicle queue has completely dissolved, as 
a function of the scaled vehicle arrival rate Qarri" for a = 
1.05. Our numerical simulation assumes the special case that 
a pedestrian arrives just at time t = to when the last vehicle 
of the queue passes the crossing point. The average time delay 
to this pedestrian is represented by triangles and compared 
to the analytical results of formula I43H (solid line). 



2. the gaps of vehicles approaching the last, already 
accelerating vehicle in a queue may be smaller than 
To. 

These two effects increase the waiting time, i.e. (ATI) < 
(A* 3 >._ 

During the waiting time (Ats + At^) of pedestrians, 
the expected number of arriving pedestrians is A (A<2 + 



At^). According to our model, all of these pedestrians 
will use the next occuring gap of size ctt or larger to 
cross the street. We can assume that the waiting time of 
the last crossing pedestrian is approximately zero, while 
it is approximately (A^2 + At^) for the first one (when 
the pedestrian arrival rate is high enough). Therefore, the 
average delay can be approximated as (At2 + Ats)/2, and 
the cumulative delay of all waiting pedestrians amounts 
to 

A ( A ^ + A ^ . (44 ) 



IV. SUMMARY AND DISCUSSION 

In this paper, we have proposed the continous-in-space 
constant-deceleration-delayed-acceleration car- following 
model (CDDA model), in order to allow for the analyt- 
ical calculation of the interactions between vehicles and 
crossing pedestrians under conditions of statistically dis- 
tributed arrival times. Although the model is not made to 
reproduce all currently known properties of traffic flows, 
it does reflect some essential features such as accident- 
free driving, a constant outflow from traffic jams and a 
characteristic queue resolution speed. 

We have distinguished two interaction modes between 
pedestrians and vehicles: (i) When pedestrians prefer 
large safety factors a > <jq, vehicles are not stopped, and 
pedestrians cross between moving vehicles either one by 
one or in small groups, (ii) When pedestrians keep small 
safety factors a < cto, they may stop vehicles, which 
usually causes vehicle queues. Once a large enough gap 
between successive pedestrian arrivals occurs, cars will 
move again and prevent the crossing of pedestrians, un- 
til the last vehicle in the queue has passed the crossing 
point. This oscillatory dynamics with alternating flows of 
cars and pedestrians tends to be inefficient and related 
with long waiting times [21). 

In this contribution, we have calculated the threshold 
(To between the oscillating and non-oscillating regime. It 
turned out to be a function of do/(ar 2 ) only, i.e. indepen- 
dent of the pedestrian or vehicle arrival rates A and Q a rr, 
the vehicle length Iq, the free vehicle velocity vo or the 
preferred time gap T, while the car deceleration a, the de- 
sired minimum distance do an d the crossing time r mat- 
ter. We have also calculated the expected waiting times of 
pedestrians and vehicles as a function of the arrival rates. 
The difficult step in gaining these results was the calcu- 
lation of the first overcritical time gap and its expected 
value. This also required the determination of the vehi- 
cle gap distribution for deterministic, i.e. non-fluctuating 
vehicle interactions, while variations in the arrival times 
were taken into account [see Eq. JJJ]. The formulas for 
the waiting time distributions can serve to judge under 
which conditions pedestrian and vehicle streams should 
be controlled (terminated) by traffic lights and when a 
self-organized crossing of streets is more efficient. Beyond 
this, our approach is generally expected to be useful for 
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a better understanding of intersecting flows and certain 
conflicting processes. For example, a similar gap accep- 
tance problem is found in lane-changing maneuvers, so 
that our formulas may help to calculate analytical for- 
mulas for lane-changing rates. 

Regarding the choice of the behavior and parame- 
ters of cars drivers (careful or aggressive) and pedes- 
trians (careful or daring), one may assume an evolu- 
tionary perspective: Due to a learning process during 
many vehicle-pedestrian interactions, an optimal behav- 
ior should emerge on the long run. It is, however, not 
yet clear whether there exists a state which is optimal 
for both, drivers and pedestrians. If not, one may con- 
sider the pedestrian-vehicle interactions as an example 
for a social dilemma |28| . and the outcome may depend 
on details of the interactions. For example, if pedestrians 
would tend to use safety factors a < o~o, car drivers may 
react to this by an aggressive approaching behavior. This 
would make it difficult for pedestrians to stop vehicles. 
However, cars could still be successfully stopped, if pedes- 
trians learned to enter a road exactly with a time gap of 



oit. In conclusion, there are always strategies to pro- 
duce or avoid alternating pedestrian and vehicle flows, 
but the outcome depends always on the parameters of 
both, pedestrian and driver behavior. The determination 
of the optimal behavioral parameters and the evaluation 
of interactive parameter adaptations of pedestrians and 
vehicles will be left for a future study. 
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APPENDIX A: CALCULATION OF THE 
EXPECTED WAITING TIME FOR A SUITABLE 
GAP 

Let P(T') be the distribution density function of vehi- 
cle time gaps T' . Moreover, let 

t, 

Q = Prob(T' < T*) = J dT' P(T') (Al) 
o 
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be the probability of finding a time gap T' < T* and 



(T') 



1 



T'<T, 



dT' T'P(T') 



(A2) 



the expected value of time gaps that are smaller than . 
Then, given that a car has just passed, the expected time 
until the first gap T' greater than T* occurs, is given by 
the expression 



(TL > ) = Y d nT'Q n {l-Q) 



(A3) 



n=0 



as an arbitrary number n of smaller gaps may occur with 
probability Q each, before a large enough gap occurs with 
probability (1 — Q). Here, we have used that the expected 
lengths T' of short gaps T' < T* just add up due to the 
assumption of independently and identically distributed 
time gaps T'. One can calculate 



{TL.) = (1-Q)T'(X>Q 

^ 17=0 



/ H 



(1 - Q)T'Q- 
QT 7 



1 



dQ VI -Q 



1-0 

J dT' T'P(T') 
o 

oo 

/ dT' p(t') 



(A4) 



Inserting the vehicle time gap distribution Q eventually 
gives 

1 - Q = P(T > T#) = (1 - g arr To)e- Qa " (T '- To) (A5) 
and 

QT ' 



Qaxr 

This implies 
(T' 



1 -(1 - QarrT )(l + Qarr T, )e" Q »' (T * " To) 

(A6) 



e Qa„(T»-T ) 
1 — QarrPo 



(1 + OarrP*) 



(A7) 



The required minimum time gap for the crossing of 
a pedestrian between two successive vehicles is T* = 
err + (Zq + do)/v, while the preferred time gap between 
successive vehicles is To = T + (Iq + do) /v. Note that 
formula l|38|l for the expected waiting time of vehicles for 
a large enough gap in the pedestrian stream corresponds 
to the special case Tq = with A = Q arr and T* = r. 



